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[57] ABSTRACT 

A method and apparatus for determining important histo- 
logical characteristics of tissue, including a determination of 
the tissue’s health. Electrical pulses are converted into 
meaningful numerical representations through the use of 
Fourier Transforms. These numerical representations are 
then used to determine important histological characteristics 
of tissue. This novel invention does not require rectification 
and thus provides for detailed information from the ultra- 
sonic scan. 

4 Claims, 5 Drawing Sheets 
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METHOD OF AND APPARATUS FOR 
HISTOLOGICAL HUMAN TISSUE 
CHARACTERIZATIONUSING ULTRASOUND 

ORIGIN OF THE INVENTION 5 

The invention described herein was jointly made by 
employees of the United States Government and an 
employee of the University of Maryland Medical School. It 
may be used by and for the Government for governmental ^ 
purposes without the payment of any royalties thereon or 
therefor. 

BACKGROUND OF THE INVENTION 

1 . Technical Field of the Invention 15 

This invention relates to the medical classification of 

human tissue as healthy or unhealthy. More specifically, the 
invention uses ultrasound for determining histological char- 
acteristics of tissue by converting the return energy pulses 
into numerical terms, thus facilitating a quantitative analysis 20 
for medical diagnosis. 

2. Discussion of the Related Art 

Ultrasound has routinely been used in the medical pro- 
fession to determine the shape, size and thickness of human 
tissue. Historically, this has been achieved through the 25 
detection of specular reflections from macroscopic tissue 
interfaces. Although this method allows for a general char- 
acterization of the tissue, it does not reveal its underlying 
health and structure. In order to obtain such a detailed 
analysis, reflections from the microscopic interfaces 3 
(scatterers) contained within the tissue must also be 
detected. An improved invention to detect these scatterers 
has yet to be discovered. 

The prior art is able to determine general characteristics of 35 
tissue from its macroscopic interface reflections. This is 
commonly determined by transforming analog ultrasonic 
radio frequency data into a visual display. This technique 
requires the expertise of a skilled operator to interpret the 
display which does not include an accurate determination of 
the underlying pathology of the tissue. 

Other prior art systems use Fourier Transforms to detect 
microscopic reflections; however, these systems do not 
allow for direct measurement of reflectance through self- 
calibration, and therefore are subject to calibration prob- 45 
lems. These methods separate the components of an ultra- 
sonic pulse into corresponding frequencies, using the 
Fourier Transform. This requires a reference plate to nor- 
malize the Fourier energy from the reflected tissue and thus 
the accuracy and range of available data is limited. 50 

Similarly, the prior art taught by Sommer, Joynt, Carroll 
and Macovski (“Ultrasound characterization of abdominal 
tissues via digital analysis of backscattered wavefronts”. 
Radiology, 141:811-7, 1981) uses a Fourier analysis to 
determine the mean spacing of scatterers in the liver and 55 
spleen. This method, however, does not provide for the 
determination of key variables that give ranges for specific 
tissues and their state of health. 

U.S. Pat No. 5,417,215 of Evans et al. provides a method 
of interpreting the microscopic interfaces which requires 60 
demodulating the return energy pulses by full-wave rectifi- 
cation to obtain the amplitude modulation of the pulses and 
computing the power spectrum by performing a Fast Fourier 
Transform on the rectified, digitized pulses. The requirement 
of rectification limits this method in providing information 65 
from the ultrasonic scan. It does not provide specific, well 
established physical parameters of the body tissue under 
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investigation. It requires the comparison of Fourier energies 
or ratios, as well as, comparisons based on correlations. 

It is, therefore, an objective of the invention to provide an 
improved ultrasound tissue characterization system for his- 
tological tissue classification. 

It is a further objective of the invention to provide a 
system which does not require rectification, thereby allow- 
ing a more detailed analysis. 

SUMMARY OF THE INVENTION 

The present invention converts electrical pulses into 
meaningful numerical representations through the use of 
Fourier Transforms. This information can be used to deter- 
mine important histological characteristics of tissue, includ- 
ing a determination of the tissue’s health. This novel inven- 
tion does not require rectification and thus provides for a 
much more detailed information from the ultrasonic scan. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is an equipment diagram for data taking; 

FIG. 2 is an illustration of an ultrasonic wave passing 
through an iterated structure constructed from two acousti- 
cally different tissue types; 

FIG. 3 is a typical ultrasonic echo pattern received from 
tissue reflections; 

FIG. 4 is a plot of the Fourier Transform equation with 
arbitrary amplitude; and 

FIG. 5 is a plot of the Fourier Transform equation with 
discrete interpolation. 

DETAILED DESCRIPTION OF THE 
INVENTION 

Referring more particularly to the drawings, a well- 
damped 10 MHz transducer one-eighth inch (Vis” ) diameter 
crystal 12 is connected to a wide band pulser-receiver 14. 
The pulser-receiver 14 is a sonic generator and receiver for 
generating, sending and receiving waves. The wide-band 
(non-rectified) output of the pulser-receiver 14 is connected 
to a digital oscilloscope 16, which has a plotter 18 connected 
to it to record the data output. The oscilloscope 16 provides 
live visual reading of the output. The data-point plotter 18 
gives a graphic reading of the output. The oscilloscope 16 
also has Fast Fourier Transform capability. 

The transducer 12 emits an ultrasonic pulse that travels 
through the body 30 into the subcutaneous tissue that is 
under examination. This tissue reflects a sequence of small 
amplitude echoes back to the transducer 12. The sequence of 
echoes from the tissue is received by the transducer 12 , 
which converts the ultrasonic pulse sequence into a 
sequence of electrical pulses. The electrical pulse sequence 
is amplified by the wide-band amplifier in the pulser- 
receiver 14 and sent to the oscilloscope 16 for measurement 
The oscilloscope 16, further processes the data, displays it 
on a visual monitor, and sends it to a plotter 18 for a 
hard-copy record. The oscilloscope 16 has the ability to 
perform Fast Fourier Transforms on the digitized pulses 20. 

Specifically, the transforming equations used are as fol- 
lows: 

1(a) 

W) = J h(l)c-Wdt = FlKt)] 
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A(f) =J' 


1(b) 


G(fyP*&dt - F[G(f)\ 


P ~ foexp 


(~S9 


cos27i/i>io 


where G(f) is the Fourier Transform of h(t). j s evaluated by substitution into equation 1(a) 

The following properties of this type of transform are 
used: 


F\h(kt)}=-j- 

here i =V-i 


exp(-i2-nft 0 )F\h(t)} (shift operation) 

2(a) 

1 n 

2(b) 

) 1 = 2 6TM01 (addition theorem) 
J m= 1 

G ^ j (similarity theorem) 

2(c) 


10 


F[ko(t)] = j Pocxp ^ -n j co&(2nf 0 f)dt 


/ > oGcxp<-7CC5 2 0 r -/o) 2 ) 

n ho(t)} 9 /> o 


( 8 ) 


(9) 
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Using standard algebraic techniques and trigonometric sub- 
stitutions the magnitude is rewritten: 


Consider the echo train from an alternating stacked sequence 
of tissues. The following relationships are used: 

P 2 ^=r<£T l2 T 2l ) 2 ^>R 12 exp(-2rtoc I / 1 _2(n-l)- 2 / 2 ) (odd mem- 
bers) 3 («) 


20 


\G(f> 0, = 0)1 : 


FoCexp(-Jio 2 (/ - /o) 2 ) 

2 


( 10 ) 


y?i2exp - 2H\ 


1- 

N T~ 


27cos(4i9fc)exi*-2/y 2 ) + (Texp(-W 2 )? 


27 2 cos(4^)exp(-2Ho) + (7 2 cxp<-2Wo)) 2 


P^~Po{T l2 T 21 f nn ' x R V 2 exp(-2n[oc 1 / 1 4« 2 / 2 ]) (even members) 3(6) 

In examining the echos, each received echo is subscripted 
with an index to indicate the temporal order of reception 

/ |0Uf] 2 \ (4) 

hj(t - to) = Pja. p / -it J 2 \ cos2 

where t 0 marks the reception of the first echo, and At j is the 
time interval between the jth received reflection and the first 
echo. 

By allowing G/f^Ffh/t-to)] and using (ltf), (2a) and 
(3tf) the following relationship is obtained 


where f 0 is the center frequency of the ultrasonic transducer, 
Hq is ajli+c^l* Hi is ctjlj, H 2 is a 2 l 2 , T is (T 12 T 21 X t 2 
is l2/c 2 , t s is 1 1 /c 1 +l 2 /c 2 , a* is the attenuation coefficient of 
the type 1 material, and a 2 is the attenuation coefficient of 
30 the type 2 material. 

An analysis of equation (10) shows the shape of IG(f>0, 
t=0)l is determined by two different undulations. The more 
rapid undulation occurs between much smaller frequency 
intervals, and is determined by the cosine function in the 
35 denominator of the radical. The cosine function, cos 4nft, 
causes this denominator to vary between the values 


G(f)- t F[hj(t-t;)) = Kmk*(t )} 

j= 1,3,5... /= 2,4,6... 

m ( ( 6 ) 

where k 1f)=Z (TnTnr'R i 2 exp ( HU + 1 F= di - V ~ 1 - 

J= 1 V 

^ (^12^21> f 'Sj^rXp ( + “2^2] “ <2 llff J — + — J ) 

It is assumed that the later reflections are small and equation 
(6) is not greatly affected by letting m approach infinity, the 
equation is rewritten: 


[l-{7* exp(-2// 0 ))] 2 , and [l+<7* exp(-2 H 2 ))f 

40 both values are perfect squares. This occurs when cosdrcfU 
varies between +1 and +1, for which dTcft^nTt or 



where Af f is the smaller frequency interval for the more 
rapid undulation. The numerator inside the radical also 
varies between 

[l-<rcxp(-2W 2 ))] 2 , and [1+<T expK^f 


* 0 ) = - 


Wuexp £ -2ai/j - i4, 9 r ” ] [ 1 ” (^i 2 ^ 2 i )exp £ -2aik - *4t/— ] ] 
[ 1 - (TuTufasp [ -21a,;, + aM +^-)] ] 


(7) 


by substituting the sum for each series. 

In general, any part of the transform may be used to 60 
express a relationship among the ultrasonic variables and, 
hence, give useful information about the significant tissue 
classifications. For purposes of illustration, the magnitude is 
used, although the real part, the imaginary part, the phase or 
any combination is just as instructive. 

Next, the Fourier transform of the incident ultrasonic 
pulse. 


but with a frequency of 



65 where M s is the frequency of the slower undulation. Solving 
equations (11) and (12) for 1 ± and 1 2 respectively it is 
determined: 
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,l= ivr 1 

' A6 1 

L 1_ TT J 

(13a) 

k= ^k 


(13b) 


By measuring the frequency intervals, Af^, and Af 2 , one can 
determine an average thickness 22 of a ground substance, 1 2 
and the cluster of substantive cells, 1 A . 

The plot of equation ( 10) reveals a series of maximum and 
minimum values. One can trace the envelopes of the relative 
maxima and minima by noting the occurrence when cos 
47tff=+l and -1, respectively 


I Guppe/j > 0, JIq — 0)1 — 


A cexpi-noHf -/o) 2 ) 
2 


x 


(14) 


constituents however it was discovered that a reasonable 
mass density model for ultrasonic applications is essentially 
a two-component mixture comprised of collagenous and 
other materials. Using the definition of mass density and 
5 further assuming that the collagen fiber diameters are 
smaller than the ultrasonic wavelength in the tissue it can be 
written 


10 


P = 


l - 




Pi 


C 

~wr 

P2 


(19) 


where C is the wet weight percent of collagen, p A is the 
density of the non-collagenous material in the tissue and p 2 
15 is the density of the collagenous material. 


/?i 2 exp(-2tfi) 


N 1 - 27exp(-2tfa)cos4TSfe + (Ttxpi-lHz)) 2 
~~ 1 - 72exp(-2/fo) _ 


while 

A aexp{-Tto 2 (f -/o) 2 ) 

I Giowdj > 0, *o = 0)1 — , y x 


20 


(15) 


N 1 - 2Texp(-2//2)cos47tffi 7 (Texpi-lfy)) 2 

/?oexp(-2//i) 

1 + T^xpi-lHo) 25 


The sound velocity in a mouse tail tendon fibers has been 
reported to be 1.733+0,056xl0 5 cm/sec (Goss and O’Brien, 
Journal of Acoustical Society of America, Vol. 65, pages 
507-511; 1979). Tendon collagen is 30% wet weight. 
O’Brien (in “Proceedings Ultrasonic International 77”, IPC 
Science and Technology Press, Guildford, England, pages 
194-205; 1977) has empirically determined that the velocity 
of sound in tissue, c. depends upon the percent wet weight 
of collagen in the tissue’s composition. C, according to the 
following relation 


For a given frequency, f 1 (e.g. at a maximum for \G uppe X^v 
t=0)l) the ratio. 


c=f 1.588+0.032 log/DxlO 5 cm/sec 


( 20 ) 


<0 = 0)1 1 + T±xp(~2H o) 

6 <o = 0)1 " i _ T^xp^^o) 

Solving for T 2 gives 

1*= 1 - «12 = ( ) e*p(2Wo) 

where e is the ratio of maximum envelope to minimum 
envelope. Using data measurements one can obtain the 
values of e as a function of frequency and by curve fitting 
obtain values for T and Hq. 40 

In determining the collagen content in the ground 
substance, equation (17) is solved for and 

Zi -7a 

Ru ~ Zi + Z 2 ’ Rn - ~ Rl1 45 


(16) 


30 


(17) 35 


is used to obtain 


Z2 =z, 




exp(2// 0 ) 


1 +\l 1 - 


e — 1 
e+l 


exp(2//o) 


(18) 

50 


where Z x is the acoustic impedance of the cluster of fat cells, 55 
which is assumed to be the same as the acoustic impedance 
of the lipids stored in the fat cell clusters. These lipids 
possess a measurable acoustic impedance. Combined with 
the experimental determination of e, Z A and Z 2 , the 
acoustic impedance of the ground substance, can be calcu- 60 
lated. This analysis suggests a closer examination of the 
factors of the ground substance’s acoustic impedance. 

In determining the collagen and mass density [gms/cm ] 
of tissues, the published values of the mass density of pure 
collagen varies from 1.16 to 1.33. However, collagen is only 65 
one of the many specialized molecular species found in 
tissues. Upon examination of the mass densities of the other 


The acoustic impedances are calculated to be 

(1.588 + 0.032 log,0 x 10 5 (21) 

7a -r; T=; gms/cm^sec 

^Tgr + jrar_ 

Pl + P2 

The measurement and analysis technique gives 26. With 
equation (21) it is possible to determine the weight percent 
of collagen in the ground substance in-vivo if the density is 
known of the collagen fibers. But the density of collagen 
fibers depends in part on the health and nutrition of the 
patient This makes the value of Z 2 a potentially good 
indicator of tissue health with the larger values of Z 2 
indicating more and better grade collagen in the ground 
substance. By assuming a reasonable range of values for the 
collagen fiber density the determination gives a range of 
values for the weight percent of collagen in tissues in-vivo 
28. The linkage is the degradation of the mechanical prop- 
erties of these tissues is directly involved with the failure of 
the tissue’s mechanical properties. That connection is the 
weight percent of collagen and the collagen’s density (the 
higher the density, the better its structural properties) in the 
ground substance. 

Using a curve fitting algorithm for the experimentally 
determined quantity 

(-*+) 

Both T 2 and Hq can be determined by solving the equation 

7l e xp(-2Ho) = (4^j-) (B) 

where H is a 1 l 1 +a 2 l 2 . Since 1 A and 1 2 are determined 
according to equations 13(a), 13(b) and a x can be estimated 
from scattering and absorption of ultrasonic waves by fat 
cells, one can determine a 2 
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H 0 - a iii (24) 

“ 2 = — V— 

which can be measured as a function of frequency. This 
information can be used to assess the quality of collagen 5 
fibers and the areal density of the fiber network (weave) in 
the ground substance. A correlation with invasive measure- 
ments and to a lesser extent to in vitro measurements will 
permit the best estimates of these quantities. 

In this manner, it is clear that more than a broad classi- 10 
fication of healthy or unhealthy can be realized for smaller 
regions of interest Most significantly, a diagnosis for the 
development of pressure ulcers can be facilitated through the 
use of this invention. 

It is specifically intended that the present invention not be 15 
specifically limited to the embodiments and illustrations 
contained herein, but embrace all such modified forms 
thereof as come within the scope of the following claims. 

What is claimed is: 

1. A method for detecting histological properties of 20 
human tissue, composed of an iterated stacking sequence of 
cellular material and collagen containing structural 
components, using data obtained from ultrasonic ray inputs, 
comprising the steps of: 

emitting a high-frequency ultrasonic pulse into human 25 
tissue; 

collecting a subsequent sequence of reflected echoes from 
said tissue; 

converting said echoes into a sequence of electrical 30 
pulses; 

amplifying said sequence of electrical pulses; 

performing a Fourier transform upon a time record of said 
sequence of electrical pulses; 

selecting any part of said transform; 

plotting a display of said transform in terms of said part; 

evaluating said display to quantitatively determine an 
index of relative collagen content in said collagen 
containing structural component, wherein said evalu- 
ating is comprised of: 
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calculating an average thickness of said cellular mate- 
rial and collagen containing components; 
determining a ratio for a series of maximum and 
minimum values that comprise said display and; 
calculating an acoustic impedance value for said col- 
lagen containing structural component. 

2. The method as in claim 1 wherein said evaluating 
includes calculating an attenuation coefficient for said col- 
lagen containing structural component. 

3. Apparatus for classifying histological properties of 
human tissue comprising: 

means for generating ultrasound pulses into human tissue; 

transducer means for applying said ultrasound pulses to 
said tissue and receiving returned ultrasound from said 
tissue body; 

means for converting said returned ultrasound pulses into 
a sequence of electrical pulses; 

means for amplifying said pulses; 

means for performing Fast Fourier Transforms on said 
amplified pulses; 

means for selecting any part of said transform; 

means for plotting a display of said transform in terms of 
said part; 

means for evaluating said display to quantitatively deter- 
mine an index of collagen content in said tissue, 
wherein said means for evaluating said display is a 
means for 

calculating an average thickness of said cellular mate- 
rial and collagen containing components; 
determining a ratio for a series of maximum and 
minimum values that comprise said display, and; 
calculating an acoustic impedance value for said col- 
lagen containing structural component. 

4. The apparatus of claim 3 wherein said means for 
evaluating said display is a means for 

calculating an attenuation coefficient for said collagen 
containing structural component. 


* * * * * 



